Josephson vortex coupled to a flux qubit 
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Experiments towards realizing a readout of superconducting qubits by using ballistic Josephson 
vortices are reported. We measured the microwave radiation induced by a fluxon moving in an 
annular Josephson junction. By coupling a flux qubit as a current dipole to the annular junction, 
we detect periodic variations of the fluxon's oscillation frequency versus magnetic flux through the 
qubit. We found that the scattering of a fluxon on a current dipole can lead to the acceleration of 
a fluxon regardless of a dipole polarity. We use the perturbation theory and numerical simulations 
of the perturbed sine-Gordon equation to analyze our results. 

PACS numbers: 74.50.+r, 84.40.Lj 
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We are experimentally investigating a new type of de- 
tector that has been theoretically proposed for very fast 
and weakly perturbing readout of superconducting qubits 
[JIH • The detection principle is based on measuring a de- 
lay time of a ballistic Josephson vortex (fluxon), moving 
in a Josephson transmission line ( JTL) [g, H[ • The vortex 
propagation delay depends on the state of the qubit mag- 
netically coupled to the line. Expected advantages of this 
method are high time resolution and weak perturbation 
of the qubit. Potentially, the method might be suitable 
for implementing quantum feedback schemes (fj @ for 
flux qubits using on-chip fluxon readout. This approach 
also offers an opportunity to incorporate the existing sin- 
gle flux quantum logic (SFQ) solutions in the readout, 
bringing the dream about scalable quantum computer 
closer to reality @, [1] ■ 

A Josephson vortex in underdamped JTL has prop- 
erties of a relativistic particle carrying a magnetic flux 
quantum $0 = h/2e [9(. The size of a vortex can vary 
from few to several hundreds of microns, depending on 
the critical current density j c and its velocity u inside 
the junction. By applying a bias current, the vortex can 
be accelerated up to the Swihart velocity cs, which is 
the speed of light in JTL. The dynamical properties of 
a fluxon resemble a classical particle with a well-defined 
mass and velocity. Nevertheless, at sufficiently low tem- 
peratures, quantum properties of fluxons such as tun- 
neling and energy level quantization have been already 
observed 

We would like to employ fluxons for developing a fast 
and sensitive magnetic field detector for measurements of 
superconducting qubits. In this Letter, we report direct 
measurements of electromagnetic radiation from a fluxon 
moving in an annular Josephson junction (A J J). The ra- 
diation is detected by using a microstrip antenna capac- 
itively coupled to the AJJ. Furthermore, we place a flux 
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FIG. 1. An annular Josephson junction with a trapped fluxon 
coupled to a flux qubit. 



qubit close to the long junction and couple them mag- 
netically with a superconducting loop (see Fig. Q]). This 
coupling scheme makes the fluxon interact with a cur- 
rent dipole [ll|, formed by the electrodes of the loop 
coupled to the qubit. The time delay of the fluxon can 
be detected as a frequency shift of the electromagnetic 
radiation emitted from the junction. This shift provides 
information about the state of the flux qubit. 

The most significant advantage of using the closed 
topology of JTL (see Fig. [2]) is the quantization of mag- 
netic flux in AJJ, which allows, by putting a pair of cur- 
rent injectors in the biasing electrode, to create a fluxon 
on demand by applying a current through the injectors 
13]. 

The circuit was fabricated using photolithography and 
standard Nb/AlO^/Nb trilayer process with the criti- 
cal current density j c ~ 1 kA/cm 2 The estimated 
Josephson penetration depth is Aj ~ 12 /mi, the Joseph- 
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FIG. 2. Optical photograph of the chip with the annular 
Josephson junction on the right part and experimental set- 
up schematics. Left part shows the zoom into the area with 
the flux qubit with a coupling loop (yellow loop) and control 
line (green loop). Red crosses indicate the positions of three 
Josephson junctions in the flux qubit loop. 



son plasma frequency uj p /2h ~ 124 GHz, and the esti- 
mated damping parameter a ~ 0.02. The circumference 
of the junction L = 1130 /im determines the frequency 
of the radiation corresponding to a single fluxon mov- 
ing with the Swihart velocity c$ to be at about 15 GHz. 
The width of the AJJ was W — 2/im and its fluxon free 
critical current I c = 23 mA. The flux qubit was made 
using the standard aluminum shadow evaporation pro- 
cess [l], UH and deposited after the niobium structures 
were fabricated. Estimated parameters for the Joseph- 
son junctions in the flux qubit loop were the following: 
critical current I c = 380 nA, alpha factor a q = 0.54, ratio 
of Josephson and charging energies Ej/Ec = 830. 

The fluxon radiation was detected first at T = 4.2 K 
temperature using of a cryogenic wide band (4-20 GHz) 
microwave amplifier with the noise temperature of about 
7 K followed by a room temperature amplifier with the 
total gain of 50 dB. The radiation spectrum was studied 
using a Rohdc&Schwarz FSUP26 spectrum analyzer. An 
example of the measured spectrum is shown in the inset 
of Fig. El Using the possibility to directly detect radiation 
of the fluxon resonant oscillations, we have performed 
systematic measurements of the dependence of the fluxon 
velocity versus bias current - the current- voltage charac- 
teristics [l6| - measured in the frequency domain (see 
Fig. [3]) . This approach provides an easy access to study 
the fine structure of the current- voltage curve as the pre- 
cision of frequency measurements is by several orders of 
magnitude greater than the resolution of direct dc volt- 
age measurement. 

To couple a flux qubit to the fluxon inside an annular 
Josephson junction, it is necessary to engineer an interac- 
tion between two orthogonal magnetic dipoles. To facil- 
itate this interaction, we have added a superconducting 
coupling loop embracing a flux qubit, as shown in Fig. [T] 
The current induced in the coupling loop attached to the 
AJJ is proportional to the persistent current in the flux 
qubit. Thus, the persistent current in the qubit manifests 
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FIG. 3. The zero-field step measured in the frequency domain 
for the ambient temperature T = 4.2 K and injection current 
Ici = 3.973 mA. The inset shows the sample spectrum for the 
fixed bias current It = 5.3 mA with the respective Lorentz fit. 



itself in the AJJ as a current dipole with an amplitude fj, 
on top of the homogeneous background of bias current. 
When fluxon scatters on a positive current dipole - it first 
gets accelerated and then decelerated by the dipole poles. 
In the ideal case of absence of damping and bias current, 
the sign of frequency change Sv is determined only by po- 
larity of the dipole. In the presence of finite damping and 
homogeneous bias current, situation completely changes - 
as the total propagation time becomes dependent on the 
complex interplay between bias current, current dipole 
strength and damping. 

A theoretical description of interaction between 
Josephson vortex and current dipole in the AJJ can be 
done by the perturbed sine-Gordon equation (PSGE) 

PCI III 



= l-M<P)+K5(*-d/2)-5(x+d/2)) 



with the periodic boundary conditions 



(1) 



Lp{-l/2,t) = <p(l/2,t) + 2irn, 



d(p(-l/2,t) _ d<p(l/2,t) 
dx 



dx 

(2) 

where n is the number of trapped fluxons, 7 = Ib/I c is 
the normalized bias current, a — u> p /uj c is the damping 
parameter, I — L/ Xj is the normalized junction circum- 
ference, fx = I^/(j c X,jW) is the amplitude of the current 
dipole and d = D/Xj is the normalized distance between 
the dipole poles. Direct analytic solution of ([l) is not 
an easy task. Therefore, we analyze Eq. ([TJ) using the 
perturbation approach developed in [l9j . In the limit of 
small perturbations 7 <C 1, ex <C 1 and motion of 

a single fluxon in the AJJ can be described by a system 
of ordinary differential equations for the fluxon's velocity 
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FIG. 4. Relative frequency deviation from equilibrium Sv/vq 
of the fluxon oscillation frequency versus bias current. Black 
line shows the result of perturbation approach, while the 
red line depicts results of direct numerical simulations of the 
PSGE equation (JSJ with = 1. The blue curve corresponds 
to the case with a<j = 0.2. 



u(t) and its spatial coordinate X(t): 



du 



^ = -^(1 - u 2 f/ 2 - au(l -u 2 )-^{\- u 2 )x 



dX 
~dt 
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(d/2 - A)sech 
(d/2 + X)sech 
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-d/2 - X 



(3) 



(4) 



By numerically solving ©-(H]) with the additional con- 
dition u(—l/2) = u{l/2) one can calculate an equilibrium 
trajectory in phase space for fluxon oscillations in the 
AJJ with the current dipole and estimate a deviation of 
fluxon oscillation frequency from the unperturbed case 



du 



vq, where vq is the oscillation frequency for 



fi = 0. Black line in Fig. @] shows the dependence of 
relative deviation Su/u versus bias current 7 calculated 
from the perturbation theory for the following set of sys- 
tem parameters: I = 20, a = 0.02, ,u = 0.05, d = 2. The 
deviation 8v is large and negative for small bias currents 
7 -C 0.1, what means that the fluxon is being slowed 
down by the current dipole and eventually can be pinned 
at the dipole if the bias current is too small. Surpris- 
ingly, for larger currents 7 > 0.05 the sign of du becomes 
positive meaning that the current dipole accelerates the 
fluxon. To understand this phenomenon, we need to look 
at the Eq. ([3]) and notice that the effective damping term 
a e = au(l — u 2 ) has a non-monotonic behavior. When 
increasing the fluxon velocity u, the effective damping 
is increasing for u < l/y/3 and then starts decreasing. 
This means that deceleration (acceleration) is favorable 
for low (high) bias currents. 



FIG. 5. Persistent current for a ground state of the flux qubit 
versus magnetic frustration (black line). Red line shows the 
corresponding fluxon shift calculated using the perturbation 
theory for. 



To verify the results of the perturbation theory, we per- 
formed direct numerical simulation of Eq. (JXJ) with delta 
functions replaced by the hyperbolic secants in order to 
smoothen current distribution: 



sech 



d 2 ip dip 

W + a ~dt 

x- d/2 



0,1 



d 2 V 
dx 2 

— sech 



7' 



sin(</?) 
-d/2 



(5) 



The parameter characterizes the width of current dis- 
tribution and is ad ~ 1 in the experiment. The prefactor 
B = l/(%ad) is chosen to keep the normalization con- 
straint B J^/^ sech((x — d/2) / a,d)dx — 1. The red curve 
in Fig.|4]shows the results of the numerical calculations of 
Eq. ([5]) with boundary conditions <j2j> - As it can be seen in 
Fig. 01 results of the numerical simulations qualitatively 
coincide with the perturbation theory. We see that this 
coincidence is improving for smaller ad as Eq. ([5]) takes 
form of Eq. ([1} in the limit of ag —> 0. 

To experimentally test the qubit readout scheme dis- 
cussed above the temperature was lowered to T ~ 70 
mK, well below the superconducting transition tempera- 
ture T c of aluminum forming the qubit. The long junc- 
tion was biased at a fixed current /{,. Then, we varied a 
current through the control line Iql m order to change 
the magnetic flux through the flux qubit. Due to the 
periodic variation of the persistent current in the qubit 
loop with the control line current IcLi the current dipole 
strength [i is modulated with an amplitude which de- 
pends on the state of the flux qubit. We can write the 
dipole strength as fi = /j,q + 5(a, where /iq is the initial 
current offset (e.g., due to trapped magnetic flux). A 
persistent current I p of the ground state in the flux qubit 
can be calculated by numerical simulation of its Hamil- 
tonian 17]. The result is depicted by the black solid line 
in Fig. E 
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FIG. 6. Modulation of the fluxon's oscillation frequency due 
to the coupling to the flux qubit. Every point consists of 100 
averages. Bias current was set at 7 = 0.521, w ~ 9.1. 



Using the presented perturbation theory (j3])-((4j) we cal- 
culated the response of the fluxon readout to the current 
dipole controlled by the signal fi = kl p (at //o = 0). The 
proportionality coefficient k is determined by mutual in- 
ductance between the qubit and the control loop as well 
as the critical current density of the AJJ. The fluxon re- 
sponse to the persistent current in the flux qubit loop 
is depicted in Fig. [5] by the red curve. It was calcu- 
lated for a fixed bias current 7 = 0.2, I = 20, a = 0.02, 
n(I p = 300 nA) = 0.05. The response is indicated in kHz 
assuming uq ~ 13 GHz (as it corresponds in experiment 
to the bias current 7 ~ 0.2). Noticeably, the response 
signal bv is approximately proportional to the amplitude 
of the persistent current and stays positive despite the 
change of the sign of I p . The asymmetry of deviation bv 
for the positive and negative branches of I p is less than 
3% . This means that the fluxon scattering is nearly inde- 
pendent of the polarity of the current dipole. It depends 
dominantly on the bias 7 and the absolute amplitude of 
the current dipole fj,. 

The experimental curve showing the reaction of the 
fluxon to the magnetic flux through the flux qubit are 



presented in Fig. [6] The periodic modulation of the 
fluxon frequency versus magnetic flux through the qubit 
corresponds to the changing of the persistent currents in 
the qubit as Fig. O suggests. We did not observe clear 
narrow peaks at the half flux quantum point, most prob- 
ably due to excess fluctuations. Emerging dip-like pe- 
culiarities can be noted at presumed half flux quantum 
points which suggest that the dips may be there, covered 
by noise and insufficient resolution. Further improve- 
ments of experimental setup are required to resolve these 
peaks. The presented measurement curve has a convex 
profile which tells that indeed the deviation of frequency 
bv is positive, consistently with predictions made above 
by the perturbation approach and numerical simulations. 

In conclusion, we have detected fluxon radiation from 
the annular Josephson junction at millikelvin tempera- 
tures. Measurements of fluxon oscillation frequency as 
a function of bias current lb resolve the fine structure 
of the current-voltage characteristics with much greater 
precision than direct voltage measurements. Using this 
technique, we have detected the modulation induced by 
the persistent current in the flux qubit coupled to the 
AJJ. We have thus implemented a microwave generator 
controlled by the flux qubit. We have observed that the 
scattering of the fluxon on a current dipole can lead to 
acceleration of the fluxon, regardless of a dipole polar- 
ity. The perturbation theory and direct numerical simu- 
lations qualitatively well describe this phenomenon. The 
tested fluxon readout scheme is compatible with Single 
Flux Quantum (SFQ) superconducting logic and can also 
be useful for applications where fast, weakly perturbing 
magnetic signal detection is needed. 
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